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ABSTRACT 

Interstellar spectra toward £ Oph acquired with the Goddard High-Resolution Spectrograph were used to 
obtain oscillator strengths for approximately two dozen S i lines. This analysis was possible because precisely 
determined experimental oscillator strengths are available for several multiplets, including one with a weak 
interstellar line. The self-consistent set of oscillator strengths then was obtained from a curve of growth based 
on line strengths spanning a range of a factor of 100. The derived /-values for a number of multiplets differ 
from values quoted by Morton (1991) but are generally consistent with the suite of available experimental and 
theoretical results. 

Subject headings: atomic data — ISM: abundances — stars: individual (£ Ophiuchi) 


l. INTRODUCTION 

The Goddard High-Resolution Spectrograph (GHRS) is 
yielding spectra of unprecedented quality. Signal-to-noise 
ratios exceeding 100-200 are routinely achieved (Cardelli & 
Ebbets 1994; Fitzpatrick & Spitzer 1994; Lambert et al. 1994). 
The resulting analyses on absorption lines are challenging the 
available atomic and molecular data from experimental and 
theoretical sources. In particular, the astronomical measure- 
ments are now requiring accuracies of 10% in more and more 
instances. The atomic data on weak intercombination lines of 
C ii, N i, O i, and Si ii (Cardelli et al. 1993, 1994) are consistent 
with the amount of interstellar absorption seen toward £ Oph 
when the comparisons include results from the damping wings 
of strong, highly saturated lines (Sofia, Cardelli, & Savage 
1994). The comparison involving strong and weak Mg n lines 
(Cardelli et al. 1991; Sofia et al. 1994), however, suggests the 
need for revised oscillator strengths for the doublet near 1240 
A. In the area of molecular /-values, several intersystem bands 
of CO detected toward £ Oph (Federman et al. 1994) indicate 
the need for revisions to the values obtained through labor- 
atory spectroscopic analyses (see Morton & Noreau 1994), and 
the observed F-X (0, 0), D-X (0, 0), and A X (2, 0) bands of C 2 
yield relative /-values that are more consistent with theoretical 
determinations than with available experimental results 
(Lambert, Sheffer, & Federman 1995). Here, data on a series of 
multiplets in S i acquired with the GHRS are analyzed so that 
a consistent set of oscillator strengths is available for future 
studies. 

Neutral sulfur has been the focus of several studies recently. 
Large-scale computations that include configuration inter- 
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action (Ho & Henry 1985; Mendoza & Zeippen 1988) have 
yielded estimates of oscillator strengths for a number of multi- 
plets. Experimental measurements have been performed on the 
relative strengths of the multiplets at 1814 A, 1479 A, and 1429 
A (Doering 1990) and on individual line /-values for the multi- 
plets >U1814, 1299 (Beideck et al. 1994). Emission from neutral 
sulfur has been reported for the Io torus (Durrance, Feldman, 
& Weaver 1983), for the vicinity of Io (Ballester et al. 1987), and 
for cometary comae (Roettger et al. 1989). Absorption from 
interstellar gas toward £ Oph has been analyzed by Federman 
et al. (1993). The astronomical data on S i are not always 
consistent with the laboratory and theoretical results, as dis- 
cussed in the above references. 

The spectral range covered by recent observations toward 
£ Oph with the GHRS span most of the multiplets for S i listed 
by Morton (1991). In this paper, 14 multiplets are analyzed in a 
self-consistent manner in order to extract a suite of oscillator 
strengths accurate to about 20% in all cases and approaching 
10% in many instances. For some, the improvement over pre- 
vious values exceeds a factor of 10. The data considered here 
include those previously discussed by Federman et al. (1993). 
The basic approach is to use the column density deduced from 
a weak line with a precise, experimentally determined /-value 
and the 6-value from analysis of the stronger lines of the multi- 
plet as a starting point for a theoretical curve-of-growth 
analysis on all the observed lines. Such an approach is possible 
because the neutral species predominantly reside in a single 
velocity component, the main component at u hclio = — 15 km 
s“\ with a well-defined 6-value (see Savage, Cardelli, & Sofia 
1992; Federman et al. 1993). This point is further highlighted in 
Figure 1, where the strong S i line at 1807 A and the Si ii line at 
1808 A are plotted on a common velocity scale. The substruc- 
ture apparent in ultra-high-resolution observations of CH and 
CN (Lambert, Sheffer, & Crane 1990; Crawford et al. 1994) 
may not be consistent with the curve-of-growth analyses for 
the atoms (Federman et al. 1993) because the two narrow 
molecular components have different strengths. From data on 
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Fig. 1. — Echelle spectra for S I X 1807 and Si n x 1 808. The strong S I line is 
only seen in the main component at — 15 km s~ l . 


K I absorption acquired at ultra-high resolution, D. Welty & 
L. Hobbs (1994, private communication) obtain results for the 
6-value consistent with those for CN and CH, but with two 
components of comparable strength. Because the K i com- 
ponents have comparable strength, a single curve of growth for 
a potentially unresolved blend should not introduce any sig- 
nificant systematic error. The reliability of our analysis also is a 
result of detections of very weak tines that lie on the linear 
portion of a theoretical curve of growth. 

2. OBSERVATIONAL DATA 

The data were acquired with the GHRS for other programs 
in which we participated. Most of the data were taken with the 
moderate-resolution grating G160M, consisting of the wave- 
length intervals 1221-1258 A, 1252-1293 A, 1296-1333 A, 
1334-1375 A, 1403-1440 A, and 1442-1488 A. Spectra in the 
vicinity of 1400 and 1807 A were obtained with the echelle 
gratings. Procedures were applied to the data to correct for 
fixed pattern noise/granularity, etc. (see Cardelli et al. 1993, 
1994; Cardelli & Ebbets 1994; Fitzpatrick & Spitzer 1994; 
Lambert et al. 1994 for details). The signal-to-noise ratios that 
were achieved in the final merged spectra (from individual FP- 
SPLITs) were between 100 and 1000. The spectra in the vicin- 
ity of each of the S I lines not already displayed by Federman et 
al. (1993) are shown in Figure 2. 

3. ANALYSIS AND RESULTS 

The standard way of analyzing data via a curve of growth 
was followed by Federman et al. (1993): weak lines set the 
column density, and stronger ones set the 6-value. The optical 
depth at line center for our data spanned a range from 0.10 to 
15. The method was applied to the data for the multiplet A 1479 
and was based on /-values recently measured in the laboratory. 
Doering (1990) obtained relative values between the multiplets 
at 1814, 1479, and 1429 A. When his ratio for/(1479)//(l 814) of 
1.12 ±0.06 is combined with the /-value for A1814 of 
0.088 + 0.005 reported by Beideck et al. (1994), a multiplet 
value of 0.099 + 0.008 is obtained. With the assumption that 
LS coupling applies, the line oscillator strengths are 
1.33 x 10“ 3 , 1.79 x 10~ 2 , and 8.03 x 10 2 for AA1474.6, 
1474.4, and 1474.0, respectively. The interstellar line at 1474.6 
A, with W x = 1.56 + 0.15 mA (Federman et al. 1993), indicates 
a column density, N ( S i), of (6.46 + 0.67) x 10 13 cm" 2 , in 
which the uncertainty in column density does not include 
uncertainties in oscillator strength. The optical depth at line 


center obtained from the analysis is 0.17, revealing that absorp- 
tion from Al 474.6 is indeed optically thin. The 6- value of 
1.1-1. 2 km s” 1 was constrained by the stronger lines in the 
multiplet. The derived 6-value is similar to the one obtained by 
Savage et al. (1992) for Mg i. [For comparison, Federman et al. 
1993 used the multiplet oscillator strength of 0.090 quoted by 
Morton 1991 for A1479, and from lines in three multiplets they 
obtained a weighted average of N( S i) = (6.95 + 0.18) x 10 13 
cm 2 and 6 = 1.2 + 0.1 kms' 1 .] 

The value for N( S i) of 6.46 x 10 13 cm -2 and a 6-value of 
1.15 km s" 1 are utilized in the present analysis for all the other 
transitions. Including the six lines (in four multiplets) studied 
by Federman et al. (1993), the present investigation is based on 
27 lines from 14 multiplets. The multiplet at 1236 A is not 
included because its lines appear on the steeply varying wing of 
a stellar feature. The results of our curve-of-growth analysis 
appear in Table 1, where the columns correspond to multiplet 
wavelength, line wavelength, the upper state designation for 
the transition (3p 4 3 P for the ground state), observed equiva- 
lent width (kFJ with 1 a measurement uncertainties, the /-value 
quoted by Morton (1991) or by Beideck et al. (1994) for the 
multiplets AA1814, 1299 (/ M ), the /-values required to place the 
observations on a single curve of growth (/ ISM ), and log/ 1SM A. 
In all cases in which the spin did not change, line oscillator 
strengths were converted into multiplet values under the 
assumption that LS coupling applies; this assumption was 
considered by Morton (1991) also and appears to be valid for 
the data obtained by Beideck et al. (1994). A visual presen- 
tation of the resulting empirical curve of growth appears in 
Figure 3. The variation about a smooth curve of growth for log 
/sm ^ — 1*5 in the bottom panel provides a measure of the 
uncertainty in 6-value from the observations. Some of the 
scatter for these values of log/, SM A arises from the estimates in 
W x for blended lines. Note that observational limitations with 
the IUE satellite prevented other empirical curve-of-growth 
analyses from including lines with optical depth at line center 
of <0.5, which is indicated in this figure. For larger optical 
depths, line saturation becomes important. 

Several points deserve mention. First, the precise experimen- 
tal /-values of Beideck et al. (1994) for lines in the multiplets at 
1814 and 1299 A also indicate 6-values of 1.1-1. 2 km s _l . 
Second, the revised oscillator strengths for the lines of the 
multiplet at 1429 A yield a multiplet /-value that agrees reason- 
ably well with the experimental results of Doering (1990) and 
Beideck et al. (1994): /, SM = 0.168 versus 0.149 + 0.011 from 
/(1429)//(1814)= 1.69 + 0.08 and /(1 814) given above. The 
line oscillator strengths for AA 1425. 19, 1425.03 are consistent 
with this assessment even though a correction to H^(1425. 19) 
of 1.90 mA was made to take into account the contribution 
from the blended weak line of the multiplet at 1425.22 A. Here, 
the equivalent width for A 1425.22 was estimated from N( S i) 
and 6 with the revised line oscillator strength from our 
analysis. This multiplet could be used in conjunction with the 
multiplet at 1479 A in future analyses if high-resolution 
(echelle) observations are available. Third, in a similar vein, the 
blended lines in the multiplet A 1320 were analyzed through 
estimates of the equivalent widths for the two weakest lines. 
Fourth, our data for the lines AA1296.17, 1295.65 had compara- 
ble errors to those reported by Morton (1975) from spectra 
obtained with the Copernicus satellite, and so our analysis was 
based on equivalent widths derived from taking weighted 
means. Fifth, slight changes in /-values for transitions studied 
by Federman et al. (1993) are the result of setting N(S i) and 6. 
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TABLE 1 


Results of Curve-of-Growth Analysis 


Am 

(A) 

X 

(A) 

Designation 

"i 

(mA) 

/m' 

f b 

J ISM 


1814 

1807.311 

4s 3 S° 

23.97 + 0.55 

9.6 X 10 2 

9.6 x 10 2 

2.239 

1479 

1474.571 

4s 3 D° 

1.56 ±0.15 

1.21 x 10" 3 

1.33 x 10" 3 

0.293 


1474.379 


11.85 + 0.20 

1.63 x 10" 2 

1.79 x 10" 2 

1.421 


1473.994 


18.78 ±0.20 

7.30 x 10 2 

8.03 x 10' 2 

2.073 

c 

1472.971 

3 d *D° 

12.88 ± 0.27 

1.91 x 10" 2 

2.18 x 10' 2 

1.507 

c 

1444.296 

4 s'D 0 

1.40 ±0.40 

8.13 x 10 4 

1.24 x 10" 3 

0.253 

1429 

1425.219 

M 'D° 

(1.90) d 

2.38 x 10 3 

1.75 x 10“ 3 

0.397 


1425.188 


14.46 ±0.31 

3.65 x 10 2 

2.69 x 10" 2 

1.584 


1425.030 


20.54 ± 0.28 

1.92 x 10" 1 

1.41 x 10" 1 

2.303 

1405 

1401.514 

5s *S° 

9.87 + 0.51 

1.61 x 10" 2 

1.49 x 10" 2 

1.320 

1320 

1316.622 

4d 3 Z)° 

(0.44) d 

4.11 x 10 4 

4.54 x 10' 4 

-0.223 


1316.615 


(5.26) d 

6.15 x 10 3 

6.79 x 10" 3 

0.951 


1316.543 


13.57 ±0.30 

3.45 x 10 2 

3.81 x 10' 2 

1.700 

1307 

1303.430 

6s i S° 

4.08 ±0.31 

2.91 x 10" 2 

5.06 x 10' 3 

0.819 

1299 

1296.174 

4s" i P v 

1 1.42 ± 0.95 e 

2.2 x 10" 2 

2.2 x 10" 2 

1.455 


1295.653 


16.99 + 0.88 e 

8.7 x 10" 2 

8.7 x 10 2 

2.052 

1274 

1270.787 

5d*D° 

(1.42) d 

1.00 x 10" 2 

1.64 x 10' 3 

0.319 


1270.780 


6.08 ± 0.21 

5.51 x 10~ 2 

9.02 x 10" 3 

1.059 


1270.769 



6.62 x 10" 4 

1.09 x 10' 4 

-0.859 

1266 

1262.860 

7s i S° 

2.15 ±0.29 


2.59 x 10" 3 

0.515 

1256 

1247.160 

6d i D° 

1.37 ±0.31 

3.24 x 10" 2 

1.64 x 10~ 3 

0.311 


1247.134 



5.77 x 10“ 3 

2.90 x 10‘ 4 

-0.442 


1247.107 



3.81 x 10" 4 

1.92 x 10“ 5 

-1.621 

1245. 

1241.905 

8s 3 S" 

0.90 ± 0.27 


1.06 x 10 3 

0.119 

1227 

1224.544 

8 d 2 D° 

<;0.80 

1.23 x 10 2 

<9.6 x 10“ 4 

<0.070 


1224.506 



2.19 x 10" 3 




1224.471 



1.48 x 10 -4 




* The results for the multiplets at 1814 and 1299 A are from Beideck et al. 1994. 
b The column density and 6-value were derived by fitting a curve of growth to the lines of the 
multiplet at 1479 A. Beideck et aPs. 1994 results for the multiplets AA 1814, 1299 are consistent with this 
curve of growth. The analysis resulted in refined /-values for the multiplets A21429, 1320, 1274, 1256 
under the assumption that LS coupling applies for line oscillator strengths. 

c Morton 1991 does not give a multiplet wavelength for the spin-changing transitions. 
d The number in parenthesis is the estimate removed from the blended feature, 
c Average of the results from Copernicus and GHRS. 


Last, the structure seen in K i absorption (D. Welty & L. 
Hobbs 1994, private communication) does not affect our con- 
clusions, in large measure because very weak lines were 
analyzed here. There is also the possibility that S i, which has a 
larger ionization potential than K l, is more widely distributed 
along the line of sight and consequently has a larger intrinsic 
b-value. The b-value of 1.5 km s _1 obtained by Federman et al. 
(1993) for the dominant ion Ni n suggests that dominant ions 
are even more widely dispersed in this picture. 

4. DISCUSSION 

As revealed in Table 1, our analysis yields significant 
revisions in /-values for 12 lines in six multiplets and two new 
determinations when comparison is made with the compilation 
of Morton (1991), updated to include the results of Beideck et 
al. (1994). These / SM -values, as well as those for the multiplets 
AA1405, 1320, and 1227, are listed in Table 2 with available 
experimental and theoretical results. Line oscillator strengths 
are only displayed when such were reported. The comparisons 
for A/U814, 1299 are discussed in Beideck et al. (1994) and are 
not repeated here. 

There is general agreement among the results shown in 
Table 2. As noted above, the experimental oscillator strengths 
of Beideck et al. (1994) form the basis for our analysis, and 
when combined with the relative /-values measured by Doering 
(1990), they provide an additional check on the value used for 


the multiplet A 1 429. Furthermore, the experimental results 
earlier than 1990 are consistent with the more accurate values 
reported here, except for Muller’s (1968) result for the line at 
1303 A. 

For the most part, there is also agreement involving the 
theoretical determinations. The correspondence between / 1SM 
and the values obtained from the large-scale computations of 
Ho & Henry (1985) and Mendoza & Zeippen (1988) is reason- 
able although not quite as good as that discussed by Beideck et 
al. The one exception is the theoretical result of Mendoza & 
Zeippen for the multiplet at 1274 A. It appears that cancel- 
lation in the dipole transition integral occurs for A 1274, a result 
reminiscent of that for /1479 (see Beideck et al. 1994), For 
/1274, the admixture of upper states involves 3p 3 ( 4 S°)5d 3 D° 
and 3p 3 ( 2 D°)3d 3 D° (Martin, Zalubas, & Musgrove 1990). The 
results from the Opacity Project (Butler, Mendoza, & Zeippen 
1995) agree with the results derived here from interstellar 
absorption remarkably well, being somewhat more consistent 
than those of Mendoza & Zeippen (1988). Discrepancies at the 
50% level still persist, however, for the multiplets AA1479, 1274, 
as well as a significant difference for X 1256. Martin et al. (1990) 
also note an admixture involving 3p\ 2 D°)3d 3 D° for ^1256. We 
remark that the results of Ganas (1982) for transitions into nd 
3 D° states differ significantly from the interstellar ones, possibly 
because configuration interaction was not taken into account. 
Moreover, the results of Kurucz and Peytremann for AA 1 274, 
1256, and 1227 quoted by Morton (1991) do not appear to be 
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Fir 3 -Curves of growth (u) based on the /-values quoted by Morton (1991) and (ft) based on the revised /-values presented here. The lines m the multiple! at 
1479 A on which the present analysis is based, are indicated by solid circles. The open circles represent lines whose oscillator strength either was based on the 
experiments of Beideck^t al. (1994) or was refined to be consistent with the curve of growth defined by .11479. The bottom panel also shows where the optical depth 
at line center is 0.5, when saturation begins to play an important role in the analysis. 


reliable; their /-values differ by factors of 5 to 20 when com- 
pared with/i SM . 

In conclusion, a curve-of-growth analysis yielded a sell- 
consistent set of oscillator strengths for 14 multiplets of S i. The 
analysis, which was based on recently obtained precise experi- 
mental values (Doering 1990; Beideck et al. 1994) for the multi- 


ple At 479, is consistent with these experimental results for 
three other multiplets, A/1814, 1429, and 1299. The four multi- 
ples are also seen in emission from solar system objects 
(Durrance et al. 1983; Ballester et al. 1987; Roettger et al. 
1989). Our analysis highlights the interplay among precise 
laboratory measurements on selected lines, high-quality 


TABLE 2 

Comparison of Results for Oscillator Strengths 


References 


ISM 

Savage & Lawrence 1966 b 

Muller 1968 b 

Doering 1990/Beideck et al. I994 b 

Lawrence 1967 c 

Gruzdev 1969 c 

Aymar 1973 (CI) C 

Ganas 1982 c 

Ho & Henry 1985 (L) c 

Ho & Henry 1985 (V) c 

Fawcett 1986 c 

Mendoza & Zeippen 1988 c 

Butler et al. 1995 c 


Reference 


ISM 

Muller 1968 b 

Smith 1978 h 

Aymar 1973 (CI) C 

Ganas 1982 c 

Ho & Henry 1985 (L) c 

Ho & Henry 1985 (V) c 

Mendoza & Zeippen 1988 c . 
Butler et al. 1995 c 


A1444 


A1425.22 


A1425.19 


A1425.03 


A1429M* 


1.24 x 10" 
8.1 x 10“ 

7.0 x 10“ 

3.0 x 10“ 
9.4 x 10" 


1.75 x 10" 


2.69 x 10“ 


2.47 


10 “ 


3.56 


1.62 


10 “ 


10 " 


1.41 x 10" 


1.94 x 10" 


9.66 x 10" 


vi 1 402 


A1405M* 


A1320M* 


A 1303 


1.49 x 10“ 
1.6 x 10" 

1.8 x 10 


1.68 x 10" 1 
(2.4 ± 0.4) x 10 1 
1.8 x 10“ l 
(1.49 ±0.11) x 10 1 


2.30 x 10" 
2.56 x 10" 
2.18 x 10" 
1.40 x 10" 


1.72 x 10" 
1.51 x 10 


A1307M* 


1.43 x 10" 2 

4.52 x 10“ 2 

5.06 x 10“ 3 

5.04 x 10 


4.2 x 10" 2 

2.9 x 10" 2 



(3.8 ±0.5) x 10“ 2 



1.8 x 10" 2 




1.5 x 10" 2 

1.15 x 10“ l 


5.4 x 10 

1.1 x 10“ 2 

3.3 x 10“ 2 



1.4 x 10“ 2 

2.6 x 10 2 



1.2 x 10" 2 

2.6 x 10“ 2 


4.0 x 10 

1.5 x 10" 2 

2.8 x 10" 2 


4.8 x 10 


Reference 

A1274M* 

A1266M* 

A1256M* 

A1245M* 

A1227M* 

ISM 

1.08 x 10" 2 

2.59 x 10" 3 

1.92 x 10“ 3 

1.06 x 10" 3 

<1.2 x 10 3 

M filler 1968 b 

7.3 x itr 3 





Ganas 1982 c 

5.68x10" 2 

2.6 x 10" 3 

3.17 x 10' 2 



Mendoza & Zeippen 1988 c 

2.0x10" 3 

2.0 x 10“ 3 




Butler et al. 1995 c 

4.3 x 10' 3 

2.1 x 10“ 3 

9.2 x 10' 5 

1.1 x 10“ 3 

1.3 x 10“ 3 


* M indicates multiplet wavelength. 

b Experimental results. _ _ .u 

c Theoretical results in which for Aymar 1973 the results based on configuration interaction are shown, and for Ho & Henry 1985 both the 

length (L) and velocity (V) results are given. 
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astronomical data for a suite of lines of varying strength, and 
large-scale theoretical computations for a large number of 
multiplets. Overall, the suite of oscillator strengths derived 
here have an accuracy of some 20%, with those for the multi- 
plets with precise experimental data approaching 10%. There- 
fore, future analyses involving lines of neutral sulfur are now 
more secure. Of particular note are those analyses of emission 
from Io and its plasma torus and cometary comae where 
uncertainties in the atomic data are no longer the limiting 
factor. 
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